Thermodynamics of Mesoscopic Vortex Systems in 1+1 Dimensions 
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The thermodynamics of a disordered planar vortex array is studied numerically using a new poly- 
nomial algorithm which circumvents slow glassy dynamics. Close to the glass transition, the anoma- 
lous vortex displacement is found to agree well with the prediction of the renormalization-group 
theory. Interesting behaviors such as the universal statistics of magnetic susceptibility variations 
are observed in both the dense and dilute regimes of this mesoscopic vortex system. 
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The behavior of vortices in dirty type-II supercon- 
ductors has been a subject of intense studies in the 
last decade [Q. Aside from the obvious technologi- 
cal significance of vortex pinning, understanding the 
physics of such interacting many-body systems in the 
presence of quenched disorder is a central theme of 
modern condensed matter physics. Similarities between 
the randomly-pinned vortex system and the more famil- 
iar mesoscopic electronic systems ^ are highlighted by 
a recent experimental study of a planar vortex array 
threaded through a thin crystal of 2H-NbSe2 by BoUe 
et al. jsj. Interesting behaviors, including the sample- 
dependent magnetic responses known as "finger prints" , 
have been observed for such a mesoscopic vortex system. 

The disordered planar vortex array is well studied the- 
oretically [|-|. It is one of the few disorder-dominated 
systems for which quantitative predictions can be made, 
including a finite-temperature "vortex glass" phase ||] 
characterized by anomalous vortex displacements Q , and 
universal variation of magnetic susceptibility How- 
ever, until the work of BoUe et al., there were hardly 
any experimental studies of this system, with difficul- 
ties stemming partly from the weak magnetic signals in 
such 2d systems. Also, numerical simulations have been 
limited by the slow glassy dynamics [0, although the 
availability of special optimization algorithms did lead 
to the elucidation of the zero-temperature problem in 
recent years In this letter, we describe numerical 

studies of the thermodynamics of the vortex glass via a 
mapping to a discrete dimer model with quenched dis- 
order. A new polynomial algorithm for the dimer prob- 
lem circumvents the glassy dynamics and enables us to 
study large systems at finite temperatures. Our results 
obtained in the dilute (single-flux-pinning) regime com- 
pare well with the experiment by BoUe et al. [Q, while 
those obtained in the collective-pinning regime strongly 
support the renormalization-group theory of the vortex 
glass, including its prediction of universal susceptibility 
variation 

The Model: The dimer model consists of all complete 
dimer coverings {D} on a square lattice £ as illustrated 



in Fig. 1(a). The partition function is 



Z = ^ exp 

{D} 



<ij>eD 



(1) 



where the sum in the exponential is over all dimers 
of a given covering, and Td is the dimer temperature. 
Quenched disorder is introduced via random bond ener- 
gies Cij, chosen independently and uniformly in the in- 
terval (-i, i). 
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FIG. 1. (a) Snapshot of a dimer covering (thick bonds) 
together with the associated "height" values h{r) for a lattice 
of size L = 8 at temperature Td — 1.0. (b) Dimer covering of 
the fixed reference, (c) Vortex line configuration (thick lines) 
obtained as the difference between (a) and (b). 

The dimer model is related to the planar vortex-line 
array via the well-known mapping to the solid-on-solid 
(SOS) model (see Figs. 1): Take the centers of the square 
cells of C to form the dual square lattice C. Orient all 
bonds of C such that the elementary squares of C that 
enclose the sites of the chosen sublattice of C (indicated 
by the solid dots in Fig. 1(a)) are circled counterclock- 
wise. It is now possible to assign a single- valued "height" 
function h{r) on the lattice points r of C, such that the 
difference of every pair of neighboring heights across the 
oriented bonds is —3 if a dimer is crossed and -1-1 other- 
wise. For the dimer covering of Fig. 1(a), the values of 
the associated height function are shown at their respec- 
tive positions. In terms of the height configuration h(r), 
the partition function (Q) can be written alternatively as 
Z = 'J2{h{r)} e~^'^^^\ where the SOS Hamiltonian takes 
the following form in the continuum limit, 



1 



(3H 



K 



{Vhf -i{r)-Vh + g cos (Gh) 



(2) 



Here, K is an effective stiffness caused by the inability 
of a tilted surface to take as much advantage of the low 
weight bonds as a flatter surface, and f(r) is a random 
local tilt bias. The periodicity of the cosine potential in 
is given by G = 27r/4 since the smallest "step" of 
this height profile is four. In the present context of a 
randomly pinned vortex array, h{r) describes the coarse- 
grained displacement field of the vortex array with re- 
spect to a uniform reference state at the same vortex 
line density; see Fig. 1(c) and Refs. 1^,^. This mapping 
between the dimer model and the planar vortex array 
allows us to compute statistical properties of the vortex 
array by monitoring appropriate quantities of the dimer 
model. The latter can be accomplished using a polyno- 
mial algorithm described below. 

The Algorithm: Computing the partition function of 
complete dimer coverings on an arbitrary weighted lattice 
is likely to be algorithmically intractable Weights 



here refer to the Boltzmann factors 



exp(-ey/Td) 



on the bonds. A weighted planar lattice Q can however 
be oriented, denote by G, so that the square root of the 
determinant of the weighted adjacency matrix defined on 
G yields the partition function (0) exactly p^ ]. Compu- 
tation of the determinant can be achieved by various row 
or column reduction schemes in time that grows polyno- 
mially with the matrix size. 

Recently, Propp and coworkers furnished the above 
algebraic reduction with an elegant graphical interpreta- 
tion for the case of bipartite planar lattices jl^. The 
partition function Zl on a lattice of linear size L is re- 
lated to Zl-i as Zl{{w}) = C{{w])Zl-i{{w']) after a 
simple weight transformation {w} — > {w'}. The pref- 
actor C{{w}) is independent of dimer coverings. The 
partition function is obtained in a "deflation" process in 
which the above recursive procedure is carried out down 
to L = with Zq = 1. This deflation process can also be 
reversed in an "inflation" process where a dimer covering 
at size L — 1 can be used to stochastically generate a dimer 
covering at size L according to Z^ already obtained. Re- 
peating the inflation process thus generates uncorrelated 
"importance samplings" of the dimer configurations, or 
equivalently the equilibrium height configurations, with- 
out the need to run the slow relaxational dynamics. The 
ensuing numerical results are obtained by taking various 
measurements of the height configurations generated this 
way. A somewhat inconvenient feature of this approach 
is that the algorithm requires open boundary condition 
on the dimer model; this in turn fixes the total number 
of vortex lines, e.g., to L/2 on a L x L lattice (see Fig. 1). 

Numerical Results: Since the temperature of the vor- 
tex array, given by K^^ in (|2|), is generally different 
from the dimer temperature Td, we first need to cal- 
ibrate the temperature scale. To do so, we exploit a 



statistical rotational symmetry |17 ^ of the system (||) 
which guarantees that at large scales, the effective K 
is the same as that of the pure system. In partic- 
ular, K can be obtained by measuring the disorder- 
averaged thermal fluctuation of the displacement field 

= "(774(o^"^^7J^^(oT2) = (27ri^)-iin(i), 

where hT{r) = h{r) — {h{r)) measures the thermal distor- 
tion superposed on the distorted "background" /i£)(r) = 
{h{r)) by the quenched disorder. We used ~< ... (...), 
and overline to denote spatial, thermal, and disorder av- 
erages respectively. With our polynomial algorithm, we 
were able to perform thorough disorder averages for equi- 
librated systems of sizes up to 512 x 512. To reduce 
boundary effects, we focus on the central L/2 x L/2 piece 
of the system and compute its displacem ent fiuct nations. 
Fig. 2(a) illustrates the dependence of W^{L) for vari- 
ous dimer temperatures Tj. The linear dependence on 
ln(L) is apparent. Identifying the proportionality con- 
stant with (27rA')^^, we obtain the empirical relation be- 
tween and Td shown in Fig. 2(b). Note that our 
result recovers the exact relation K^^jTd oo) = 16 /tt 
for the dimer model without disorder pq | . Since the glass 
transition of the system is expected to occur at tem- 
perature Kg'^ — 47r/G^ = 16/7r, our system is glassy for 
the entire range of dimer temperatures. 
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FIG. 2. (a) The disorder-averaged thermal displacement 
fluctuation W^{L) for various dimer temperatures Td. Each 
data point was obtained by generating 10^ thermal samplings 
for each of 10'^ disorder realizations. The straight lines are 
fit to W| = {2TTKy^ ln(L) + const, (b) The extracted rela- 
tion between and Td. The horizontal line indicates the 
asymptotic value of A'~^ as Td oo. 



A striking feature of the vortex glass phase is the 
anomalous fluctuation of the disorder-induced displace- 
ment, W^) = (-< hjj{r) > < /i_D(r) ^'^). The renormal- 

ization group (RG) theory predicts that W^|,(L) ~ 
C2ln^(L) for large L, with the proportionality constant 
C2 depending quadratically on the reduced temperature 
r = Kg/ K — 1 just below the glass-transition temper- 
ature K~'^ . We show W^{L) vs. ln(L) for different 

temperatures in Fig. 3(a). The forms of W^{L) are 
well described by a quadratic function of ln(L) as pre- 
dicted. The solid lines are least-square fit to the form 
wEiL) = C2 ln2(L) + Ci ln(L) + Cq. 
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FIG. 3. (a) The anomalous displacement fluctuation Wl, 
at various dimer temperatures averaged over 10'' disorder re- 
alizations. The lines are fits to the quadratic dependence on 
ln(Z/) expected from the RG theory, (b) Log-log plot of the 
quadratic coefficient C2 (extracted from (a)) vs. the reduced 
temperature r = Kg/K — 1. The best power law fit yields an 
exponent of 2.04(9). 

To probe the RG prediction further, we examine the 
temperature dependence of the quadratic coefficient C2. 
Plotting C2 extracted from the quadratic fit in Fig. 3(a) 
vs. the reduced temperature t = Kg/K — 1 using the 
vortex temperature K^^ obtained in Fig. 2(b) and the 
exact glass-transition temperature K^^ = IS/tt, we ob- 
tain the data shown in Fig. 3(b). A quadratic depen- 
dence of the coefficient C2 on the reduced temperature 
T is clearly demonstrated by the data, thereby providing 
strong support for the RG theory of the vortex glass. 

Another interesting feature of the vortex glass phase is 
the sample-to-sample variation of the magnetic suscepti- 
bility X- It was predicted |^ that the fractional variance, 

var[x]/x' for < -r « 1, (3) 

is universal, with D being a computable, size- 
independent constant of order unity. Due to the con- 
straint of fixed vortex density when using the dimer 
representation, it is not easy to probe the magnetic 
susceptibility directly by varying vortex density. In- 
stead, we use the fluctuation-dissipation relation, e.g., 
X oc {{dxh)^) — {dxh)'^ . To circumvent lattice effects, 
we take as definition x(L) = q^((/iq/i_q) — (^q)(^-q)), 
where /iq is the Fourier transform of h{r), and the right- 
hand side is evaluated at {q^, Qy) = (1/i, 0). As shown in 
Fig. 4, the fractional variance is indeed size-independent 
over the range L = 32 to L = 256 examined. Its tempera- 
ture dependence can again be deduced from the relation 
between K~^ and Td given in Fig. 2(b); it is well de- 
scribed by the linear form (^) close to the glass transition 
with D = 0.454(5). 

The universal variation of magnetic susceptibility is 
reminiscent of the phenomena of universal conductance 
fluctuation much studied in mesoscopic electronic sys- 
tems. This mesoscopic vortex system in fact exhibits 
interesting fluctuations even deep into the dilute regime, 
i.e., close to the threshold field Hd- But to study the 
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FIG. 4. Fractional variance as a function of the reduced 
temperature, r, for L = 32, 64, 128, and 256 averaged over 
10^ disorder realizations. The solid line is a least-square lin- 
ear fit for |r| < 0.15. 

vortex behavior there, we must relax the constraint of 
fixed vortex number. To do so, we focus on a Lx x Ly 
located at the center of the L x L lattice. Random bond 
energies are assigned only to this "inner" region. While 
the total number of vortex lines must be fixed, the num- 
ber of lines in the inner region can be changed by adding 
an extra bond energy Seij = {—iy~^^v in a staggered fash- 
ion to all vertical bonds of the staggered reference dimer 
pattern (Fig. 1(b)) in this region. This has the effect of 
increasing the "self-energy" of the vortex lines such that 
when V is sufficiently large, i.e., for v > Vc, all vortex lines 
will be expelled from the inner region to the outer region; 
see Fig. 5. A dilute regime in the inner region is then 
achieved for small Sv = Vc — v. This setup could be inter- 
preted as embedding the dirty type-II superconductor of 
focus within another pure type-II superconductor which 
has a much smaller Hd- In the limit L ^ L^, the outer 
region can be regarded as an infinite reservoir in which 
the vortex density change is negligible. The response we 
measure will thus come essentially from the inner region. 




v=0.00 v=0.32 v=0.40 

FIG. 5. Snapshots of vortex line configurations in the inner 
region at different repulsion strength v for a given disorder re- 
alization and Td = 0.01. Quenched disorder exists only in the 
inner region of size Lx = Ly = L/2 = 64. 

In Fig. 6(a), we plot the thermal averaged number of 
vortex lines (n) in a thin strip of size Lx — 16, Ly — 256 
embedded in a 512 x 512 square lattice. We vary the 
repulsion strength u, holding the temperature fixed at 
Td — 0.01 (corresponding to only ^ 10% of the glass tran- 
sition temperature). The solid squares are data taken 
from a given realization of disorder. We see a disordered 
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stair-case structure in the dependence of (n) on dv, sim- 
ilar to experimental magnetic response curves observed 
by Bolle |^ . The step width dispersion results from 
a delicate energetic balance between the line-line repul- 
sion and the random pinning; it can be regarded as a 
finger print of this disordered sample. By averaging such 
stair-case-response curves over 100 realizations, we ob- 
tain results indicated by the open squares in Fig. 6(a). 
We can extract a linear dependence, i.e., (n) ^ Sv (ex- 
cept for finite-size rounding very close to the threshold). 
This finding agrees well with the experimental result of 
Bolle et al^. It should be contrasted with the very dif- 
ferent behavior when the disorder is absent. The latter 
is shown in Fig. 6(b) for a system which is twice larger in 
all linear dimensions. The data are well described by a 
square-root singularity, (n) ^ (Sv)^^'^, as expected from 
the theory of Pokrovsky and Talap ov [T^ . The origin of 
the anomalous linear dependence of (n) on 6v in the ran- 
dom case has been discussed previously it can be 
readily understood from the properties of a single vortex 
line pinned in a two-dimensional random potential. 




(a) " (b) 

FIG. 6. (a) Average number of vortex lines (n) vs. difler- 
ent repulsion strength v at temperature Td — 0.01 in a narrow 
Lx X Ly stripe. The solid squares are for a given disorder real- 
ization and the open squares denote results from average over 
100 realizations. The line is a fit to linear dependence in the 
vicinity of the threshold Vc- Slight rounding near Vc is due 
to the small but still finite ratio L^/ L. (b) in) vs. v in the 
absence of disorder. The line is a fit to the square-root form 
expected in the vicinity of Vc- 

Conclusion. Finite temperature simulations were 
performed on a disordered dimcr model to study vari- 
ous equilibrium properties of a planar array of vortex 
lines. Numerical results on the anomalous displacement 
fiuctuation strongly support the renormalization-group 
theory of vortex glass. Universal susceptibility variations 
in the collective pinning regime were observed in accor- 
dance with the theory; critical behaviors near the glass 
transition compared well with analytic predictions. Suit- 
able modification of the dimcr model allowed a direct 
study of the dilute regime where a few vortex lines pen- 
etrate the system. Disordered stair-case-like magnetic 
finger prints were obtained; the statistics are in agree- 
ment with theoretical predictions, and drastically differ- 
ent from the well-known behavior of the pure system. 



The numerical findings in the dilute regime are also in 
agreement with recent experiment on a micron-sized sin- 
gle crystal of 2H-NbSe2. It is hoped that the present 
study will stimulate further experimental investigations 
of the fascinating physics of mesoscopic vortex systems. 
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